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ABSTRACT 

We combine a deep 0.5 deg^, 1.4 GHz deep radio survey in the Lockman Hole with infrared and 
optical data in the same field, including the SERVS and UKIDSS near-infrared surveys, to make the 
largest study to date of the host galaxies of radio sources with typical radio flux densities ^ 50 /rJy. 
87% (1274/1467) of radio sources have identifications in SERVS to AB « 23.1 at 3.6 or 4.5/xm, and 
9% are blended with bright objects (mostly stars), leaving only 4% (59 objects) which are too faint 
to confidently identify in the near-infrared. We are able to estimate photometric redshifts for 68% 
of the radio sources. We use mid-infrared diagnostics to show that the source population consists of 
a mixture of star forming galaxies, rapidly accreting (cold mode) AGN and low accretion rate, hot 
mode AGN, with neither AGN nor starforming galaxies clearly dominating. We see the breakdown in 
the K — z relation in faint radio source samples, and show that it is due to radio source populations 
becoming dominated by sources with radio luminosities ~ 10^^ WHz“^. At these luminosities, both 
the star forming galaxies and the cold mode AGN have hosts with stellar luminosities about a factor of 
two lower than those of hot mode AGN, which continue to reside in only the most massive hosts. We 
show that out to at least 2^2, galaxies with stellar masses > 10^^'^ Mq have a radio-loud fraction 
up to ~ 30%. This is consistent with there being a sufficient number of radio sources that radio-mode 
feedback could play a role in galaxy evolution. 

Subject headings: galaxies: active ~ radio continuum: galaxies - galaxies:evolution 


1. INTRODUCTION 

The radio source population at 1.4 GHz flux densities 
5'i.4 < 1 mjy contains a mixture of active galactic nuclei 
(AGNs), both radio-loud and radio-quiet, and star form¬ 
ing galaxies (Gondon 1992, Afonso et al. 2005, Simpson 
et al. 2006). Radio-loud AGNs have, historically, been 
easier to study, being three to four orders of magnitude 
more luminous in the radio than starforming galaxies and 
dominating radio surveys at bright flux density levels. 
However, recent studies have shown that within popula¬ 
tions with characteristic flux densities of tens to hundreds 
of microJanskys (hereafter the pJy population), the pre¬ 
dominance of radio-loud AGNs gives way and the popula¬ 
tion begins to include significant numbers of radio-quiet 
AGNs and starforming galaxies (Jarvis & Rawlings 2004; 
Simpson et al. 2006; Huynh et al. 2008; Smolcic et al. 
2008; Seymour et al. 2008; Strazzullo et al. 2010; Mao et 
al. 2012; White et al. 2012; Bonzini et al. 2013). Ghar- 
acterizing this mixed population is important for using 
radio surveys to study the star formation history of the 
Universe (Karim et al. 2011; Zwart et a. 2014), the cos¬ 
mic history of AGN activity (Smolcic et al. 2015), and 
in order to use radio continuum surveys for cosmological 
purposes. 

Multiwavelength observations of radio sources allow 
the study of the relative distributions and properties of 
the three dominant types of sources (radio-loud AGN, 
radio-quiet AGN, and starforming galaxies) among the 
various wavelengths. Previous studies have used this 
multiwavelength technique to examine radio source pop¬ 
ulations in several deep fields, but have reached differ¬ 
ing results as to the contributions of the different galaxy 
types. Most studies suggest a significant fraction (Huynh 
et al. 2008; Simpson et al. 2006; Smolcic et al. 2008) or a 


majority (White et al. 2012; Seymour et al. 2008) of star¬ 
forming and radio-quiet galaxies within the /rJy popula¬ 
tion, although some authors argue a continuation of the 
dominance of radio-loud AGNs from the Jansky and mil- 
liJansky populations to sub-milliJansky levels (Mignano 
et al. 2008). 

Deep radio surveys, combined with multiwavelength 
data, are essential for studying the nature of the host 
galaxies of the AGN-powered radio source population, 
and, as they are deep enough to contain a mixture of both 
radio-loud and radio-quiet objects, are helpful in our un¬ 
derstanding of the origin of radio-loudness in AGN. At 
high 1.4 GHz radio luminosities (L1.4 > 10^^ WHz“^) 
radio sources have long been known to be hosted by 
the most massive galaxies at all redshifts. Early work 
by Lilly & Longair (1984) and Lilly (1989) showed that 
the near-infrared iL-band magnitudes of radio galaxies 
formed a tight Hubble diagram. This so-called “K — z 
relation” has been found to be a very useful tool for esti¬ 
mating the redshifts of distant, luminous radio sources, 
and is only very weakly dependent on radio luminosity 
for radio sources with L1.4 > 10^^ Wm“^ (Willott et al. 
2003; McLure et al. 2004). The scatter in the AT — z re¬ 
lation towards fainter AT-magnitudes has long been pre¬ 
dicted to increase as fainter radio samples are studied, 
however (e.g. Jarvis et al. 2001), when radio source pop¬ 
ulations become dominated by progressively lower radio¬ 
luminosity AGN and starforming galaxies. The first signs 
of this breakdown in the K — z relation start to appear 
in p-Jy radio surveys (Simpson et al. 2012), which we 
confirm with a larger sample of sources in this paper. 

Understanding the mix of galaxies that host the /iJy 
radio source population is essential for planning for fu¬ 
ture radio observatories such as the SKA, in particu- 
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lar for cosmological studies where different types of host 
galaxies cluster differently (e.g. Ferramacho et al. 2014). 
We therefore decided to use some of the best currently- 
available multiwavelength data to investigate the nature 
of the radio source host population. We base our mul¬ 
tiwavelength catalog on the Spitzer Extragalactic Rep¬ 
resentative Volume Survey (SERVS), which uses deep 
Spitzer data to detect galaxies out to z ~ 5 (Mauduit 
et al. 2012). All five fields chosen by this survey were se¬ 
lected specifically to be integrated into multi-wavelength 
observations. In particular, the Lockman Hole field has 
two deep radio surveys within it, the very deep survey of 
Owen et al. (2009) in the north of the field and the wider, 
but somewhat shallower survey of Ibar et al. (2009) in the 
east. In this paper, we match the SERVS data with the 
Ibar et al. radio survey, and other surveys in the Lock- 
man Hole field at optical through infrared wavelengths 
to investigate the host galaxy properties of the ^Jy ra¬ 
dio source population. These data are summarized in 
Section 2. Details of the catalogs and band-matching 
are given in Section 3, and Section 4 discusses the qual¬ 
ity of the photometric redshifts. The methodology for 
discriminating between AGN and starforming galaxies is 
presented in Section 4. Section 6 presents the properties 
of the radio sources, and conclusions may be found in 
Section 7. A cosmological model of Dm = 0.3, Da = 0.7, 
and Hq = 70 kms“^Mpc“^ is assumed throughout the 
text. 

2. MULTIWAVELENGTH OBSERVATIONS 

The Lockman Hole is a particularly appealing candi¬ 
date for a multiwavelength study because there have been 
many observations made of it in a wide variety of wave¬ 
lengths. These data can be combined to form a very 
thorough study encompassing objects with redshifts out 
as far as z Ri 5. We have assembled a multiwavength 
dataset, detailed below, with which to study the host 
galaxies of /iJy radio sources in this field. Figure [1] shows 
the spatial overlap of the datasets described below. 

2 .1. SERVS data 

SERVS uses deep Spitzer data taken during the 
postcryogenic Spitzer Space Telescope mission over the 
course of 1400 hours using the Infrared Array Camera 
(IRAC) at wavelengths of 3.6 and 4.5 ^m, and detects 
objects out to z « 5. The survey covers a total of 18 deg^ 
to a detection limit of ss 2 iiJy [AB = 23.1) in five fields, 
including 4 deg^ in the Lockman Hole. SERVS catalogs 
were produced using Sextractor (Bertin & Arnouts 1996) 
as detailed in Mauduit et al. (2012). 

2 .2. Optical data 

Gonzalez-Solares et al. (2011) observed the Lockman 
Hole with the Isaac Newton Telescope Wide Field Cam¬ 
era (INT WFC) in the u,g,r,i, and z-bands, with AB 
magnitude limits of 23.9, 24.5, 24.0, 23.3, and 22.0, re¬ 
spectively. Supplementary observations were taken with 
the Mosaic 1 camera on the Mayall 4-m Telescope of 
the Kitt Peak National Observatory (KPNO) in the g, r, 
and i-bands. Details of the source extraction are given 
in Gonzalez-Solares et al. (2011). 

2.3. Ground-based near-infrared data 
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Fig. 1. — Surveys in the Lockman Hole used in this work. The 
greyscale shows the SWIRE 24/i.m image, which is close to coinci¬ 
dent with both the IRAC SWIRE coverage, and the optical imaging 
survey of Gonzales-Solares et al. (2011). The black outline shows 
the coverage of the UKIDSS DXS survey, the blue the SERVS cov¬ 
erage and the red circles the positions of the VLA pointings used 
for the Ibar et al. (2009) radio survey. The green square shows 
the deep (“level 3”) Herschel HerMES coverage (the whole field is 
covered to a lower depth), see Section 5.5. 

The UKIRT Infrared Deep Sky Survey (UKIDSS) ob¬ 
served the Lockman Hole as part of its Deep Extra- 
galactic Survey (DXS) using the Wide Field Camera 
(WFCAM) of the United Kingdom Infrared Telescope 
(UKIRT) at wavelengths 1.2 fim and 2.2 p,m, corre¬ 
sponding to the J and K bands, respectively. The sur¬ 
vey reaches AB magnitude depths of Jab = 23.1 and 
Kab = 22.5. Lawrence et al. (2007) describe the survey 
and the production of the catalogs. 

2.4. SWIRE data 

The Spitzer Wide-area Infrared Extragalactic 
(SWIRE) survey observed the Lockman Hole dur¬ 
ing its cryogenic mission using the Infrared Array 
Camera (IRAC) and the Multiband Imaging Photome¬ 
ter for Spitzer (MIPS). IRAC observed the Lockman 
Hole using all four channels, 3.6 p,m, 4.5 p,m, 5.8 ^m, 
and 8.0 p,m; and MIPS observed the Lockman Hole with 
all three of its channels, 24 ^m, 70 /rm, and 160 /rm. The 
SWIRE catalogs reach depths of ~ 10 /rJy {AB « 21.4) 
at 3.6 and 4.5 ^m. SWIRE detected galaxies out to 
z Ri 3 (Lonsdale et al. 2003). Details of the data 
release used may be found in the SWIRE data release 2 
documenlQ. 

2.5. Radio data 

Part of the Lockman Hole was observed at frequen¬ 
cies of 610 MHz and 1.4 GHz using the Giant Metre- 
wave Radio Telescope (GMRT) and the Very Large Ar¬ 
ray (VLA), respectively by Ibar et al. (2009). Their 
survey covered k. 0.5 deg^ in the Southern region of 
the Lockman Hole (centered on RA, Dec 163.0°,57.35°, 
generally called the “Lockman Hole East”) and contains 
1452 radio sources detected at 1.4GHz, comprised of 1467 
discrete components. The GMRT data reach an RMS 
noise of « 15 /rJy/beam, and the VLA data an RMS of 

^ http://irsa.ipac.caltech.edu/data/SPITZER/SWIRE/docs/delivery_doc_r2. 
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Ri 6 /ijy/beam in the center of the surveys, where the 
sensitivity is highest, tapering off to 10 /xJy/beam at 
1.4 GHz towards the edge of the survey region (see Fig¬ 
ure [2] for the flux density distribution).The depth of the 
GMRT data is such that nearly all of the 1.4GHz sources 
have counterparts at 610MHz. 

3. CATALOG MATCHING 

Data from the optical through far-infrared surveys 

were combined into a SERVS “Data Fusion” catalog, 

" 

using a matching radius of 1.0 to match them to 
the SERVS dual-band catalog. This catalog contained 
951102 objects covering the full 4 deg^ of the SERVS 
Lockman area. Details of this catalog will be presented 
in Vaccari et al. (in preparation). This Eusion catalog 
was then matched to the 1467 individual component de¬ 
tections at 1.4 GHz (with peak-to-noise ratio > 5) in 
the Ibar et al. (2009) catalog. The radio sources were 
matched to the SERVS Fusion catalog using a selec¬ 
tion radius of 1.2 arcseconds (compared to a 4 arcsecond 
beamsize), which resulted in 1245 matched objects - a 
85% match rate down to the AB Ri 23.1 catalog limit. We 
then compared the remaining 222 unmatched radio com¬ 
ponents with the SERVS data manually, and matched an 
additional 29 radio components, six of which we believe 
are components of multi-component radio sources, and 
23 of which were matched to faint SERVS sources below 
the catalog limit, resulting in a 87% match rate. Of the 
remaining unmatched objects, 134 (9% ) were confused 
by brighter objects (mostly stars and bright galaxies) in 
SERVS. Of the remainder, only 31 (2%) appeared to be 
completely blank fields, and a further 28 were most likely 
blank, but were within « 2" of faint objects in the field. 
These statistics emphasise the rarity of infrared-faint ra¬ 
dio sources (Norris et al. 2011). 

To evaluate the number of matches expected at ran¬ 
dom, we offset the positions of the radio sources by ap¬ 
proximately 1-arcmin in RA and Dec, and redid the 
matching four times with different offsets. The mean 
number of random matches was 90, corresponding to 6% 
of the number of radio components, which we treat as 
negligible. 

4. PHOTOMETRIC REDSHIFTS 
4.1. Determination of Photometric Redshifts 

Photometric redshifts of the objects observed in the 
SERVS study were obtained by Pforr et al. (in prepa¬ 
ration) using the Hyper-z photometric redshift code 
(Bonzella, Miralles & Pello 2000). A correlation of the 
photometric redshifts with the spectroscopic redshifts of 
AGN in the field (see below) is described in the Ap¬ 
pendix, and shows that ^ 85% of the redshifts are ac¬ 
curate to within an RMS of Az/(1 -|- z) Ri 0.06, the re¬ 
mainder being outliers. Highly-luminous AGN at high 
redshifts, where the hot dust from the AGN in the IRAG 
3.6 and 4.5/im bands can affect the redshift estimate, 
may have poorer photo-zs, however. 

Of the photometric redshifts for the radio galaxies, 20 
(1%) failed to find a solution. All of these failed solutions 
had only limits on their detection in either or both of the 
J-band or /F-band. To ensure a good quality photomet¬ 
ric redshift, we thus excluded the 254 objects (20% of 


the 1245 objects matched between the radio and SERVS) 
lacking a J- or iA-band detection from any analysis that 
needed redshift information, ensuring that we had a min¬ 
imum of four bands to estimate a redshift (J, 7A, [3.6] 
and [4.5]). In Figure|2|we show the 1.4 GHz flux density 
distribution of the survey for objects with and without 
photometric redshifts, indicating that there is not a very 
strong bias towards, for example, only the brighter radio 
sources having photometric redshift information. The 
photometric redshifts extend out to z « 5, as seen in Eig- 
ure|3| (though some of the redshifts at z > 4 use only the 
minimum four bands, and so may be unreliable). There 
is a dip in the redshift distribution at z ss 1.4 which has 
no clear origin. Two possible explanations are that it 
may be a degeneracy in the photometric redshifts in the 
region of the “redshift desert” where the optical bands 
in particular do not contain strong spectral features to 
constrain the redshifts, or that it may reflect real large- 
scale structure in the field. Spectroscopic redshifts were 
available for 62 of the AGN from Lacy et al. (2013) or 
the SDSS quasar survey (Ahn et al. 2014), including the 
normal type-1 AGN whose optical emission is dominated 
by quasar light. We have used these in preference to the 
photometric redshifts where applicable. 

4.2. Predictions of Photometric Redshifts for ^Jy radio 

sources 

We used the SKA Simulated Skies (S^) extragalactic 
model (Wilman et al. 2008) that uses measured lumi¬ 
nosity functions combined with an underlying numerical 
dark matter simulation to obtain a predicted redshift dis¬ 
tribution for fiJy radio sources. This predicted redshift 
distribution is compared to our measured one in Figure 
O Both distributions peak at z ~ 0.7, though our mea¬ 
sured distribution has fewer sources in the lowest redshift 
bin. This deficit is also seen in the smaller (but much 
more spectroscopically-complete) sample of Simpson et 
al. (2012), so is unlikely to be due to large-scale structure 
fluctuations, though we cannot rule out incompleteness 
at very low redshifts due to galaxies being resolved out 
in the radio survey. We see fewer sources in the z > 2.5 
bins than the predictions, this may be explained though 
by the 20% of faint sources we excluded due to miss¬ 
ing J- and/or AT-band detections. The median redshift 
we observe, 0.99, is very similar to the median redshift 
predicted by the simulations (1.04). 

5. OBJECT CLASSIFICATION 

We used a combination of a mid-infrared color-color di¬ 
agram and the radio to mid-infrared flux ratio to classify 
as many of the radio sources as possible into three overall 
classes, namely cold-mode AGN (high accretion rate ob¬ 
jects with dynamically-cold accretion flows and Edding¬ 
ton Luminosities pEdd ^ 0.01 — 1), hot-mode AGN (low 
accretion rate objects accreting gas at about the galaxy 
virial temperature of ~ 10®K, with LEdd 5 0.01, and 
jets possibly powered by the spin of the black hole, e.g. 
Martinez-Sansigre & Rawlings [2011]) and star-forming 
galaxies. Mao et al. (2012) show that mid-infrared clas¬ 
sifications are broadly consistent with those from spec- 
trosopy. 


^ http://mattiavaccari.net/df 


5.1. Rapidly accreting, cold-mode AGN 
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When looking for the relative contributions of AGN 
and starforming galaxies, bright AGN are simpler to 
isolate than starforming galaxies. For the 479 ob¬ 
jects that were detected in all four of the 3.6, 4.5, 
5.8 and 8.0/rm bands of the Infrared Array Gamera 
(IRAG) in the SWIRE data, we were able to isolate 
rapidly-accreting, cold-mode AGN using the method pro¬ 
posed by Lacy et al. (2004, 2007), which selects AGN 
based on their color relative to starforming or quies¬ 
cent galaxies. Cold-mode AGN have a characteristic red 
power-law in the mid-infrared that stands out well in a 
logio(S'5^S'3.6) versus logio(5'8.Q/5'4,5) color-color plot, 
(Figure lH left). We used Sss and ,54.5 from SERVS 
and /Fs.s and /S'g.o from SWIRE. Three cuts were made: 
logio(5'5.8/<S'3.6) > —0.1, logio(<S'8.o/<S'4.5) > —0.2, and 
logio(5'8.o/-S'4.5) < 0.81ogio(5'5.8/53.6), which produce 
the wedge-shaped area in the upper right corner of Fig¬ 
ured] This area represents the objects likely to be cold¬ 
mode AGN. These 226 objects make up 15% of the radio 
sources in the Ibar et al. survey (FigurejS]). Spectroscopic 
follow-up of AGN samples selected in this manner show 
that this particular “wedge” selection is reliable in the 
sense that 78% (527 out of 672) of objects for which spec¬ 
troscopic redshifts and classifications could be obtained 
using optical or near-infrared spectra are confirmed as 
AGN (Lacy et al. 2013). (Many of the remainder may 
also be highly-obscured AGN, but cannot be confirmed 
as such from optical/near-infrared data alone.) Among 
the cold mode AGN there are a total of six type-1 quasars 
from the Sloan Digital Sky Survey Data Release 10 (Ahn 
et al. 2014) (five of which are in the Lacy et al. [2013] 
sample). 

5.2. Low accretion rate, hot-mode radio AGN 

Lacy et al. (2006), Afonso et al. (2011), Gurkan et 
al. (2014), and Singh et al. (2014) also find that some 
weak radio AGN appear outside the AGN or starforming 
galaxy selection region, with IRAG colors similar to qui¬ 
escent galaxies. These galaxies are likely the low accre¬ 
tion rate, “hot mode” AGN, whose radio jets are highly 
efficient. These radio sources can be found in the area be¬ 
neath the cut formed by the line logio(58,0/54.5) = —0.2 
in the left-hand panel of Figure H) 

We also use the 924 = logio(524/5i4oo) value to classify 
objects without 4-band IRAG detections (e.g. Ibar et al. 
2008). As the SWIRE 24/rm flux density limit, ~ 150/rJy, 
is five times 5i.4 ghz of our faintest radio sources, and 
given that the typical value of <724 for star forming galax¬ 
ies is ss 1 (Appleton et al. 2004; Ibar et al. 2008), we ex¬ 
pect to detect the majority of star-forming galaxies and 
cold-mode AGN at 24 /rm. We indeed find that the frac¬ 
tion of radio sources detected at 24/rm is high (overall, 
874/1245 (70%) of matched radio sources are detected at 
24 ^m). Thus, the radio-bright objects not detected at 
24/rm are likely to be hot-mode AGN. We place the di¬ 
viding line between radio-loud and radio-quiet at (724 ~ 0 
at 2 = 0, then use the template fc-corrections shown in 
figure 1 of Ibar et al. (2008) to approximate the change 
in redshift of this dividing line as (724 < — 0.15(1 -I- z) (the 
dashed line in Figure 1, right-hand panel). For those ob¬ 
jects not detected at 24 /rm we assume a limit of 150 T^Jy 
to classify them. We separate the hot mode AGN se¬ 
lected using the wedge-based and (724-based methods in 
Figures |3| and [H but thereafter include them all in a 


single class, that comprise about 4% of the total radio 
source population in this sample. 

For the objects lacking photometric redshifts we ap¬ 
proximated a classification by assuming a mean redshift 
of unity (dotted line in the right-hand panel of Figure]!]). 
As can be seen, this should result in < 10 misclassihca- 
tions. 

5.3. Star-forming galaxies and 24 fJ-m detected objects 
with ambiguous classifications 

The remainder of the objects in the infrared color-color 
diagram - the upper left corner in Figure]!]- correspond 
to starforming galaxies, and make up 15% of the radio 
sources. In addition, 30% of objects are detected at 247im 
with (724 above the dashed (or dotted, in the case of ob¬ 
jects lacking photometric redshifts) line in the right-hand 
panel of Figure ]!] but not detected in all four IRAG 
bands. These most likely comprise a mix of cold-mode 
radio AGN and star-forming galaxies, and we henceforth 
classify them as “ambiguous AGN/SF”. 

5.4. Unclassifiable objects 

32% percent of the identified radio sources remain un¬ 
classifiable using either the infrared color-color diagram 
or (724, because they lack 4-band IRAG detections, are un¬ 
detected at 24 fim and their flux density limit at 24 fxm 
is insufficiently deep compared to their radio fluxes to 
classify them as radio-loud or radio-quiet. 

5.5. Comparison of wedge and (724 selection 

For the most part our classification scheme is self- 
consistent between the wedge technique and the (724 one. 
Amongst the objects with spectroscopic redshifts, how¬ 
ever, 15/37 of the hot mode AGN that are classified as 
such by the wedge technique show up as radio-quiet in 
the (724 plot. Furthermore, in eight cases, objects initially 
classified as hot-mode AGN are detected in Herschel Her- 
MES (Oliver et al. 2012) images at 250 /im and also have 
photometric redshifts > 1. These are most likely mis- 
classified star-forming Ultraluminous Infrared Galaxies 
(ULIRGs), which, due to their high redshifts, no longer 
have the PAH features in the IRAG 87rm band resulting 
in them falling out of the star formation region in Eig- 
ure]!] (left). They have been reclassified as star-forming 
galaxies. There is also one low redshift (z = 0.04) galaxy 
which appears to be a HerMES detection, we also reclas¬ 
sify it as a starburst. The remainder tend to plot in the 
radio-louder side of the infrared-radio correlation. They 
may have significant residual star formation, or 24:fj,m 
emission from their AGN (although hot-mode radio AGN 
tend to be faint in the mid-infrared, some do have signif¬ 
icant mid-infrared emission, e.g. Ogle, Whysong & An- 
tonucci [2005]). Similarly, there are 4/207 galaxies clas¬ 
sified as star forming in the wedge selection that appear 
radio-loud in the (724 one, these may be radio-loud AGN 
in starburst galaxies. Further, more accurate, classifica¬ 
tion could be attempted with full SED fitting, including 
Herschel data from the HerMES survey (Oliver et al. 
2012), but we will defer this activity to a future paper. 

Our classification results can be compared to the study 
of Bonzini et al. (2013), who used a survey of similar 
depth in the Extended Ghandra Deep Eield South to 
classify a sample of 883 radio sources as 19% radio-loud 
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Fig. 2.— Histogram of flux densities at 1.4 GHz for our sample. The dark grey indicates objects with > 4-band photometric redshifts, 
the lighter grey objects lacking photometric redshifts. 




Fig. 3.— Left: photometric redshift distribution by type (see Section 5). Blue indicates actively accreting (cold mode) AGN; dark magenta, 
hot-mode AGN (based on the IRAC color-color plot); light magenta, likely hot-mode AGN (based on 924 ); red, ambiguous AGN/starforming 
objects, and cyan, unclassifiable objects. Right: the red histogram shows the normalized photometric redshift distribution for our sample, 
the cyan one the predicted redshift distribution from the full volume of the SKA simulations (Wilman et al. 2008). 
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Fig. 4. — Classification of radio sources. Left: infrared color-color diagram. The four lines separate the objects into three areas: the 
wedge-shaped area in the upper right corner contains likely cold-mode AGN (blue dots); the area occupying the lowest two corners contains 
likely hot-mode AGN (black dots); and the area in the upper left corner contains likely starforming galaxies (yellow dots). Note there are 
a few objects classified as starbursts in the hot-mode AGN region, these are Herschel detections and are mostly z > 1 ULIRGs (see Section 
5.5). Right: 924 = logio('S'24/5'i4oo) (not fc-corrected) against redshift for the radio sources with matches at 24/im in the SWIRE catalog. 
Downward pointing triangles indicate upper limits. The dashed line indicates the criteria we used to separate radio-loud hot-mode AGN 
from cold-mode AGN and starbursts for objects with photometric redshifts, the dotted line for those without (see text for details). 
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Fig. 5. — Pie chart for the whole Ibar et al. 1.4 GHz survey using our classifications described in Section 5. The numbers of objects in 
each class are as follows; 226 cold mode AGN, including six type-1 quasars; 74 hot mode AGN (comprising of 37 from the 4-band IRAG 
classification and 37 from the ^24 classification); 207 star-forming galaxies; 437 objects ambiguous between star formation and cold-mode 
AGN; 464 unclassifiable objects (excluding those that are definitely unidentified in SERVS), and 59 unidentified objects. 


AGN, 24% radio-quiet AGN and 57% starforming galax¬ 
ies. The main difference between our classification tech¬ 
niques and those used by Bonzini et al. is their use of 
a less conservative 924 criterion for defining radio-loud 
AGN (called “hot mode” AGN in this paper), and also 
their X-ray data and deeper Spitzer data that allows 
them to resolve more ambiguities between cold-mode 
AGN and starforming galaxies. 

6 . PROPERTIES OE RADIO SOURCES 
6.1. Radio luminosity and spectral index 

In Figure El we plot the radio luminosity versus the 
photometric redshifts of the radio sources. The ra¬ 
dio luminosity was calculated for objects within the 
^Jy population using the flux densities and spectral 
indices of the 1.4 GHz objects within the Ibar et al. 
(2009) survey. Most of our radio sources have L1.4 ^ 
® WHz“^, spanning the traditional bound¬ 
ary in luminosity between radio-loud and radio-quiet ob¬ 
jects. With few exceptions (probably mis-classifications) 
all the star-forming galaxies selected via Figure 1 have 
Ti ,4 < WHz“^, corresponding to star formation 

rates ;$ 2000 assuming a Salpeter IMF. 

The Ibar et al. catalog includes spectral indices be¬ 
tween 610MHz and 1400MHz (af^Qp) for objects detected 
at both frequencies. The mean a^^QQ for the objects 
matched to SERVS is -0.78. The mean spectral indices 
for the cold-mode AGN and star forming galaxies are 
very similar to the overall mean (-0.77 and -0.79, re¬ 
spectively). That for the hot-mode AGN is signihcantly 
flatter, -0.44, indicating a large fraction of flat-spectrum 
AGN in that class. 

6 .2. K-band emission 


AT-band emission from radio galaxies tends to consist 
primarily of stellar light (Lacy et al. 1995; Best, Lon- 
gair & Rottgering 1998; Simpson, Rawlings & Lacy 1999; 
Seymour et al. 2007), although in very luminous cold¬ 
mode AGN there can be significant contributions from 
reddened quasar continuum (Rawlings et al. 1995) or 
emission lines (Bales & Rawlings 1993). Previous results 
generally found that the K — z relation corresponded to a 
curve representing a stellar population which formed its 
stars at very high redshifts {z > 5) and passively evolved 
afterwards (e.g. Lacy, Bunker & Ridgway 2000; Jarvis et 
al. 2001; Willott et al. 2003). In Figure[7]we plot the K- 
magnitudes against the logarithms of their photometric 
redshifts. The dashed curve in Figured whose function 
is K = 17.37 -I- 4.53 logio(2;) — 0.31 (logio(2;))^ , was pro¬ 
posed by Willott et al. (2003) as the polynomial of best 
fit for luminous steep-spectrum radio galaxies within a 
standard 63.9kpc metric aperture. Willott et al. show 
that their K — z relation fit closest to a model in which 
most stars formed in a single burst at redshifts of z = 10, 
and evolved passively afterwards. However, the scatter 
towards brighter magnitudes at z > 2 apparent in both 
Willott et al. and in Figure [3 cannot eliminate a model in 
which many radio source hosts formed at a lower redshift, 
z ~ 5. This is discussed further in Section 6.3 below. 

Figure [3 shows that most of the hosts of ^Jy radio 
sources are fainter in A'-band than the mean for the 
hosts of the much radio-brighter sources from the 3G, 
6 G and 7G samples (with typical radio flux densities at 
1.4 GHz ^ 50 — 1000 mJy) studied by Willott et al. 
(2003). This is perhaps not surprising for the starform¬ 
ing galaxy population, which has a very different phys¬ 
ical mechanism for radio emission than the AGN pop¬ 
ulations that dominate the Willott et al. sample, but 
it also seems true for the cold-mode AGN population. 
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at least at low redshifts/luminosities. In Figure [8] we 
plot the offset from the K — z relation for different types 
of object (especially AGN) as a function of radio lumi¬ 
nosity and redshift. In the lowest radio luminosity bin, 
2 ^q 22.5-23.5 the mean and error in the mean of 

the offset is 0.93±0.06 for cold-mode AGN (significant at 
15 (t), but insignificant (0.05 ± 0.11) for hot-mode AGN. 
In the next bin, 1023-5-24.5 the offsets is barely 

significant (0.22±0.07 for the cold-mode AGN and again 
insignficant 0.10±0.16 for the hot-mode AGN), and simi¬ 
larly we see no significant offsets in the highest luminosity 
bin (1023-5-24.5 wHz-i). 

Figure [9] shows the Petrosian magnitudes 0 of the ra¬ 
dio sources overlaid on the K — z distribution from the 
entirety of the SERVS Fusion catalog, including optical 

3 Petrosian magnitudes include the light out to a radius where 
the surface brightness of the galaxy falls to 0.2 of its mean value 
within that aperture. 


and infrared sources which were not matched to radio 
sources. Despite the offset from the K — z relation for 
the most luminous radio galaxies that we observe, the 
radio matched objects are, at any given redshift, among 
the most luminous objects, and are generally much more 
luminous than the objects found only in the optical, in¬ 
frared or ultraviolet surveys. 

6.3. Stellar mass 

We estimated the approximate stellar masses of the 
Lockman Hole radio sources using the 3.6 fim flux density 
and a 1.4 Gyr stellar population template from Bruzual 
& Chariot (2003) to calculate the rest-frame AT-band lu¬ 
minosity, Lk- The 3.6/im wavelength is a good compro¬ 
mise between being close to rest-frame AT-band at z ^ 1, 
while at the same time not usually containing significant 
contamination from hot dust emission from the AGN 
(Caputi 2013). The choice of stellar population template 
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Fig. 7.— The K — z relation, split by radio source type as defined in Section 5. The dashed line indicates the K — z relation for radio 
sources from much brighter surveys as defined by Willott et al. (2003). The hot-mode AGN are the only class of object that appears to 
trace the same K — z relation as samples of radio sources selected from surveys ~ 1000 times brighter in radio flux density. 


makes only a small difference to the /c-corrections in the 
near-infrared. We then used iC-band mass-to-light ra¬ 
tio estimates at 2 = 0 — 1 from Drory et al. (2004) and 
at z 2 from Borys et al. (2005) to parameterize the 
evolution of the iC-band mass-to-light ratio in a typical 
massive galaxy as M/Lk ~ 1-0/(1 + 1.5z). The result is 
shown in Figure[Tni Some of the largest indicated masses 
for the cold-mode accretion AGN are probably spurious, 
produced by contamination of the 3.6^m flux density by 
AGN hot dust emission. For the remainder, we see a 
fairly well-defined cutoff at 10^^'^ Mq. 

We further estimated the stellar masses of the opti¬ 
cal and infrared sources which were not matched to ra¬ 
dio sources, plotting the distribution of stellar masses 
of the radio-loud (Li .4GHz > 10^^ WHz ^) and radio¬ 
quiet populations in Figure [5] This Figure shows that 
the radio-loud sources tend to be among the objects with 
the largest stellar masses at any redshift see also Bronzini 
et al. (2013), who find that the radio-loud AGN in their 


work are also in host galaxies significantly more massive 
than cold-mode AGN or star forming galaxies, though 
the distribution of masses of our “hot mode” AGN lacks 
the tail to low masses (below ss Mq seen in figure 

9 of Bronzini et ah. A high mean mass for the radio-loud 
population is expected, based on the work of Best et al. 
(2005) and Simpson et al. (2012). 

Effective AGN feedback depends on a high duty cy¬ 
cle of activity (Groton et al. 2006), so a high fraction 
of radio-loud AGN is needed if radio-jet powered AGN 
feedback models are to be plausible. Among objects with 
stellar mass greater than or equal to 10^^ ® Mq with red- 
shifts 0.5 < z < 2.5, we found 19/69, or 28%, of all mas¬ 
sive galaxies are radio-loud (Figure [TT]), consistent with 
this requirement. At 2; > 2.5 the space density of galax¬ 
ies with masses IQ^^'^Mq is very low, ^ 10“® Mpc“^ 
(Duncan et al. 2014), and remains similar to the space 
density of radio sources. It is thus quite plausible that 
when a galaxy (and its black hole) reaches a certain crit- 
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Fig. 8 . — Offsets from the K — z relation of Willott et al. (2003) as a function of redshift and radio luminosity (symbol colors as in 
Figure 4). Top panels, all radio source types. Bottom panels, AGN-powered radio sources only. The large grey and cyan points represent 
averages, with 3 a error bars, of the hot-mode and cold-mode AGN, respectively. 


ical mass it is able to switch on as a radio-loud object 
with a duty cycle of a few tens of percent. Better statis¬ 
tics are needed at high-z, however, both on the radio 
source population and on the numbers of massive galax¬ 
ies to understand the evolution of radio source hosts at 
the earliest epochs. 

The star-forming galaxies also have relatively high stel¬ 
lar masses, ~ 10^^ Mq. The likely explanation for this is 
that these star-forming galaxies are at the high end of the 
radio luminosity function for star-forming objects, with 
typical radio luminosities ^ WHz“^, forming stars 

at rates 10MQyr“^ (Chabrier IMF; Karim et al. 2011). 
Karim et al. (2011) and Zwart et al. (2014) show that, at 
the redshifts at which these objects are mostly seen in our 
study (z 0.3 —0.5), these objects have a typical specific 
star formation rate SSFR ~ 0.1 Gyr“^. Although less 
massive objects tend to have higher SSFR, the trend 
is relatively weak, so less massive galaxies would be ex¬ 
pected to have lower absolute star formation rates and 


thus be missing from the survey. 

7. GONCLUSIONS 

The results presented in this paper include one of the 
first estimates of the redshift distribution and AT — z re¬ 
lation for a survey with a large fraction of radio sources 
< 100 ^Jy. The redshift distribution we find features 
a median value of z Ri 1 and a tail out to redshifts of 
z w 5. This is not entirely consistent with the predicted 
redshift distribution of ^Jy radio sources from the S^ 
simulations, which suggests a larger fraction of sources 
at z < 0.25, although the median redshift is very similar, 
and issues of over-resolution in the radio may explain the 
apparent deficit of very low-z sources. 

The galaxies detected in the deep radio survey at 
/ijy levels tend to be more luminous in iF-band than 
objects which were not detected in the radio survey. 
Nevertheless, we see evidence that the K — z relation 
is breaking down, not only due to the increased frac- 
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Fig. 9.— The K — z relation by type for the radio sources, with 
the overall A'-magnitude distribution for the SERVS-Lockman data 
fusion catalog shown as a log-scaled density field. Symbols as in 
Figure 3. The K — z relation from Willott et al. (2003) is shown 
as the dashed line. 

tion of star-forming galaxies, but also because the cold 
mode AGN radio sources with 1.4 GHz radio luminosi¬ 
ties Li .4 ^ 10^“^ WHz“^ tend to exist in less luminous 
host galaxies. Thus lower luminosity radio sources are 
no longer found in only the most massive galaxies, and 
so, as Simpson et al. (2012) point out, the AT — z re¬ 
lation defined for powerful radio sources should not be 
used for estimating redshifts in ^Jy surveys. It is possi¬ 
ble that the L1.4 ;$ 10^'’’ WHz“^ cold-mode AGN have a 
different mechanism for producing radio emission rather 
than AGN jets. Star formation (Kimball et al. 2011), or 
synchrotron emission powered by shocks in thermal out¬ 
flows (Zakamska & Greene 2014) have been suggested 


as possible alternate mechanisms that dominate over jet 
powered radio emission at lower luminosities. Hot mode 
AGN, on the other hand, continue to be found in only 
the most massive galaxies at Li,4 < 10^'*’ WHz“^ and 
most likely continue to be powered by jets. This picture 
is consistent with only the most massive galaxies contain¬ 
ing black holes that are capable of producing the most 
powerful radio jets, and thus with a strong dependence 
of jet luminosity on black hole mass amongst AGN (e.g. 
Laor 2000; Lacy et al. 2001; Kratzer & Richards 2014). 
We show that the fraction of massive galaxies that are 
radio-loud continues to be high (up to ~ 30%) out to 
z ^ 2, a prerequisite for models where feedback from 
radio jets is able to limit star formation in their host 
galaxies, extending the results of Best et al. (2005) and 
Simpson et al. (2012). 
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APPENDIX 

ACCURACY OF THE PHOTOMETRIC REDSHIFTS 

Spectroscopic redshifts for some of the AGN in this field were obtained by Lacy et al. (2013). We selected all 53 
of the type-2, red type-1 (defined in Lacy et al. 2013 as having broad emission lines, but optical colors redder than 
normal quasars) and “non-AGN” objects (i.e. AGN candidates selected in the mid-infrared, but which showed no AGN 
signatures in their optical spectra) in common to the two samples, with spectroscopic redshift quality g = 1 in Lacy et 
al. (2013). The plot of spectroscopic versus photometric redshift is shown in Figure[l2l The two sets of redshifts agree 
well, though the increasing dominance of the AGN-component in the higher luminosity AGN above 2; « 1 introduces 
more scatter at high redshifts, and may also be responsible for some of the catastropic outliers (15% of the matched 
objects). The scatter A2;/(l + z) Ki 0.06, though there is a small offset for the type-2 objects of Zphot — ^^spec ~ 0.15 
around z = 0.4, probably due to emission line contamination of broad-band magnitudes. The type-2 objects in this 
plot are among the most infrared-luminous of the “cold-mode” AGN in this paper, so will be the most prone to 
contamination by AGN flux in the near-infrared. 







